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ABSTRACT 

In this paper, we report our studies on the gaseous and chemical properties of a relatively large 
sample (53 members) of blue compact dwarf galaxies (BCDs). The results of correlations among the 
oxygen abundance, stellar mass, gas mass, baryonic mass, and gas fraction are present both for E- and 
I-type BCDs, which are classified according to Loose & Thuan (1985) and show elliptical and irregular 
outer haloes, respectively. These correlations of I-type BCDs show similar slopes to those of E-type 
ones. However, in general, E-type BCDs are more gas-poor and metal-rich than I-type ones at a given 
baryonic mass. Based on these results, we suggest that E-type BCDs, at least a part of them, and 
I-type ones might be likely at different evolutionary phases and/or having different progenitors. Our 
investigation of the correlation between oxygen abundance and gas fraction shows that BCDs appear 
to have not evolved as isolated systems, but to have experienced some gas flows and/or mergers. 
Subject headings: galaxies: dwarf - galaxies: abundances - galaxies: evolution 



1. INTRODUCTION 

Dwarf galaxies are by far the most numerous galaxies 
in the universe (e.g. Mateo 1998) and play a crucial role 
in our understanding of the formation and evolution of 
galaxies. They are proposed to be the building-blocks in 
hierarchical scenarios of galaxy formation (Kauffmann 
et al. 1993). Dwarf galaxies are typically metal-poor 
(e.g. Mateo 1998), and have low stellar masses. The low 
metallicities of these systems make them having simi- 
lar conditions of high-redshift galaxies. Therefore, dwarf 
galaxies at the local universe can be used as a proxy 
to reveal the physical properties of some galaxies at the 
early universe. 

Based on their morphological properties and lumi- 
nosity, dwarf galaxies can be classified into different 
types, e.g., blue compact dwarfs (BCDs), dwarf ellipti- 
cals (dEs), dwarf irregulars (dls), and etc. BCDs, which 
are the focus of the current paper, have dramatically 
different properties compared to normal dwarf galaxies 
(Zwicky 1966; Gil dc Paz ct al. 2003, hereafter G03). 
They are small, gas-rich (Hi mass fraction typically 
higher than 30%; please see Section 2.3 and 3.1) extra- 
galactic objects. Being metal-poor (0.03Z© < Z < Zq, 
here 12 + (O/H) = 8.69 (Allende Prieto et al. 2001; 
Asplund et al. 2009); e.g. Terlevich et al. 1991; Izotov 
et al. 1994; Izotov & Thuan 1999; Hunter & Hoffman 
1999), these galaxies are allowing us to study the star 
formation and chemical enrichment in a near-primordial 
environment. As such, they are important to our under- 
standing of galaxy formation and evolution. 

BCDs are spectroscopically characterized by a faint, 
blue optical continuum accompanied, in most cases, by 
intense and narrow emission lines, due to the intense 
star formation activity distributed in one to several star- 
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forming regions (Cairos et al. 2001a). dEs have smooth 
visual appearances and low gas content, and they repre- 
sent the bright, gas-poor dwarf galaxies which have little 
star formation (e.g. Binggeli 1994). By contrast, dls arc 
characterized by a lack of spiral arms and an irregular 
visual appearance (see review by Hunter 1997). They 
have low dust content (e.g., Calzetti 2001) and large gas 
reservoirs (see, e.g., Hoffman ct al. 1996). Central sur- 
face brightness increases with increasing luminosity in 
all of these types of dwarf galaxies (see recent review of 
Tolstoy et al. 2009), thus it seems natural to consider 
a possible evolutionary relationship between BCDs, dls 
and dEs (e.g. Thuan 1985; Ferguson & Binggeli 1994; 
Papaderos et al. 1996). 

However, the relation between BCDs, dls, and 
dEs/dSphs remains unclear. Current knowledge bearing 
on this question comes mainly from analysis of the pho- 
tometric properties (such as structure, surface brightness 
profile, and color) of these three types of dwarf galaxies 
(e.g. Thuan 1985; James 1994; Papaderos et al. 1996; 
Noeske et al. 2000; Gil de Paz & Madore 2005; Vadu- 
vescu et al. 2006). The results from different authors 
seem to be contradictory. For example, Thuan (1985) 
suggests that BCDs and dls are the same type of galaxy, 
differing only in their present star formation rate, and 
they can evolve into dEs. However, James (1994) obtains 
an opposite result that BCDs are probably dEs under- 
going starburst activities, while dls are a fundamentally 
distinct population from them. The results in Papaderos 
et al. (1996) show that BCDs have a brighter central 
surface brightness and a smaller exponential scale length 
than dls and dEs at a given B-band luminosity. However, 
Gil de Paz & Madore (2005) claim that at least -15% of 
BCDs in their sample have structural properties compat- 
ible with being dEs (also see Papaderos et al. 1996) ex- 
periencing a burst of star formation, whereas Vaduvescu 
et al. (2006) conclude that BCDs have similar structures 
to dls. Recently, Vaduvescu & McCall (2008) suggest 
that there are strong connections among BCDs, dls, and 
dEs since these galaxies lie precisely on a fundamental 
plane defined by the K s -hand absolute magnitude, the 
central surface brightness, and the neutral hydrogen line 
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width (also see McCall et al. 2012). 

Some other attempts to address this issue rely on the 
chemical evolution of dwarf galaxies. Richer & McCall 
(1995) show that [O/Fe] for a number of field dEs are 
larger than those in field dls at similar luminosities and 
conclude that there is unlikely an evolutionary connec- 
tion between dEs and dls. Vaduvescu et al. (2007; here- 
after VMR07) argue that the evolution of BCDs has been 
similar to that of dls since they have the same metallicity- 
gas fraction correlations. However, the sample in VMR07 
is too small (14 BCDs), and therefore might suffer large 
uncertainties. Moreover, the correlation between metal- 
licity and the gas fraction for BCDs and dls shows large 
dispersions (see Figure 10 in VMR07). 

In our previous works, we have studied the luminosity 
(or mass)-mctallicity (L-Z) relation (Zhao et al. 2010; 
hereafter ZGG10), stellar population and star formation 
(Zhao et al. 2011; hereafter ZGG11) for a sample of 
BCDs. Our results show that BCDs are old systems cur- 
rently experiencing starbursts, and follow a well-defined 
L-Z relation. Based on these two works and G03, along 
with Hi data compiled from literature (see reference in 
Table 1), we can further investigate the evolution of 
gaseous and chemical properties of BCDs. In this paper, 
we will explore correlations among different galactic pa- 
rameters, i.e. oxygen abundance, stellar mass, gas mass, 
baryonic mass, and gas mass fraction. We also com- 
pare these relations for morphologically different types 
of BCDs and dls. Our investigations on such correla- 
tions may shed some light on the answer to the question 
of whether there is an evolutionary connection between 
BCDs, dls and dEs. 

The paper is organized as following: Section 2 de- 
scribes the sample of BCDs and the data. In Section 
3, we first calculate the gas masses using the compiled 
Hi fluxes and derive the stellar masses, and then present 
several relationships involving stellar mass, gas mass and 
oxygen abundance. We also do some comparisons be- 
tween dls and BCDs in this section. Based on the results 
shown in Section 3, we discuss possible evolutionary con- 
nections between different types of dwarf galaxies, and 
gas flows in BCDs, in Section 4. A brief summary is 
given in the last section. 

2. SAMPLE AND DATA 

2.1. The Sample 

Using a unified concept of BCD galaxy, G03 have com- 
piled a BCD sample (including 105 members) from sev- 
eral exploratory studies. The galaxies in their sample 
were selected by putting forward a new set of quantita- 
tive classification criteria. This new set of criteria is a 
combination of the galaxy's color, morphology and lumi- 
nosity. Briefly, a BCD galaxy has to fulfill the following 
three criteria (referring to blue, compact and dwarf, re- 
spectively): (1) ^s.poak - Mi?,poak < 1, where fi B , pca k 
and /U_R lP eak are the peak surface brightness of B- and R- 
band, respectively, (2) ^B,pcnk < 22 mag arcsec -2 , and 
(3) the absolute magnitude Mk > —21 mag (see G03 for 
details) . 

To study the chemical evolution of BCDs, we need to 
know their metallicities and gaseous properties. Based 
on the sample in G03, ZGG10 have selected a subsam- 
ple of BCDs whose metallicities have been measured or 



having SDSS spectroscopic data. Therefore, the sam- 
ple we used here is derived from ZGG10 and a galaxy 
is selected if its Hi data has been given in literature. 
The final sample contains 53 BCD galaxies, whose basic 
properties are listed in Table 1. The distance, B and R 
photometric data and morphological type (except three 
galaxies; see below) are all adopted from G03. Except 
for several galaxies whose distance had been determined 
by measuring the magnitude of the tip of the red giant 
branch (RGB), the majority were computed using the 
galactic standard of rest velocity of the galaxies assuming 
a Hubble constant, Hq, of 70 km s _1 Mpc" 1 (Frccdman 
et al. 2001). For the photometric data, the foreground 
Galactic extinctions have been corrected in G03, using 
Ab values determined following Schlcgcl et al. (1998) 
and the Galactic extinction law of Cardelli et al. (1989) 
with Ry = 3.1. We did not attempt to correct for in- 
ternal extinction, since this is known to vary spatially 
in BCDs (e.g., Cannon et al. 2002; Garcfa-Lorenzo et 
al. 2008; Cairos et al. 2010). The internal extinction 
obtained through Balmcr lines and/or continuum fitting 
based on the aperture spectra can only be applied to a 
small part of a galaxy. In addition, it was derived only 
for about 40% of the total sample. 

The morphological types of these BCDs were classi- 
fied according to Loose & Thuan (1986). However, for 
Mrk 005 (Cairos et al. 2001b; Amorin et al. 2009), Mrk 
178 (Schulte-Ladbeck et al. 2000; Noeske et al. 2003) 
and Mrk 1329 (Noeske et al. 2003; Vaduvescu et al. 
2006), we adopt their morphological types as iE, based 
on the deeper optical and/or IR observations in the above 
references. Please refer to G03 for details of the afore- 
mentioned properties. In the following, we refer to all 
nE and iE BCDs (showing an outer elliptical envelope) 
collectively as "E-type" BCDs, and refer to all il BCDs 
(showing an irregular outer halo; here including two iO 
BCDs) as il or "I-type" BCDs unless otherwise stated. 

The oxygen abundances are taken from ZGG10 (also 
see references therein). In the majority of the cases the 
oxygen abundances have been determined on the basis of 
electronic temperatures measured from the [O in] A4363 
line (i.e. T e -method). In a few cases T e -based abun- 
dances are not available, we adopt the N2-based oxygen 
abundances, which have been derived with the empiri- 
cal N2-0/H relation proposed by Denicolo et al. (2002) 
(i.e. N2-method, where N2= log ([Nil] A6584/H Q )). As 
shown in ZGG10, for the BCD sample studied here, 
metallicities derived with the N2-method are generally 
consistent with those derived with the T e method, and 
have an accuracy of ~ 0.2 dex (also see Denicolo et al. 
2002; Pettini & Pagel 2004 and Salzer et al. 2005). 
Please sec ZGG10 and references therein for details. 

In Figure 1 we plot the integrated (B — R) color and the 
absolute B magnitude (Mb) distributions for the G03 
sample and our subsamples using different histograms. 
We can see that ours have similar distributions of (B — R) 
and Mb to the G03 sample, except that our il sub-sample 
contains a higher fraction of galaxies with Mb ~ 16.5—17 
mag and a lower fraction of galaxies with (B — R) > 1.0. 
Therefore, our subsamples can be representative of the 
total sample of G03. However, please note that none of 
these samples are complete, although they covers a large 
range of galactic parameters. 

In order to explore whether their is an evolutionary 
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TABLE 1 
Properties of our sample BCDs 
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References. — (1): Gordon &: Gottesman 1981; (2): Thuan & Martin 1981; (3): Huchtmeier & Richter 1988; (4): Thuan et al. 1999; (5): 
Salzor ot al. 2002; (6): Gavazzi et al. 2005; (7): Huchtmeier et al. 2005. 

Note. — Column 1: galaxy name; Column 2: distance moduli. Column 3—4: B andi?,-band apparent magnitudes, respectively. Column 5: 
morphological type. Column 6: oxygen abundances adopted from ZGG10; the error for the N2-bascd mctallicity is the typical scatter comparing 
to those obtained with T e method. Column 7: method of determining oxygen abundance. Column 8: total gas mass. Column 9: stellar mass. 
Column 10: method of determining stellar mass (please sec Section 3.1 for details). Column 11: gas fraction. Column 12: inverse gas fraction. 
Column 13: references for Hi data. 

a All data, except for morphological types of three galaxies (Mrk 005, Mrk 178 and Mrk 1329), arc adopted from Gil de Paz et al. (2003), and 
all of the magnitudes have been corrected for the Galactic extinctions. 
^The stellar mass is derived with: 1— i^ s -band data; 2— integrated spectra; 3— SDSS fiber spectra; 4— i?-band data. 
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Fig. 1. — Distributions of absolute B magnitude (upper) and 
(B — R) color (bottom) for the three BCD samples. Our subsamplcs 
(hatched filled and solid open) show similar distributions to the 
G03 sample (dashed open). 

link between BCDs and dls by means of comparing their 
gaseous and chemical properties, we also include a sam- 
ple of field dls adopted from Lee et al. (2003) in our 
analysis. Lee et al. (2003) sample is compiled from the 
literature and consists of 22 dls. The data for this dl 
sample, including metallicity, gas mass and stellar mass, 
are all adopted directly from Lee et al. (2003). The 
mctallicities are determined by the direct method, and 
the gas masses are derived using the same method as 
described below. Whereas for the stellar mass, we will 
check the consistency in the next section. 

2.2. Stellar Mass 

The stellar mass, M*, of a galaxy, can be derived from 
its luminosity and mass-to-light ratio (M*/L), which can 
be estimated by comparing observed galaxy colors (e.g. 
Bell & de Jong 2001), broadband spectral energy distri- 
butions (e.g. Walcher et al. 2008), or full optical spectra 
(e.g. Cid Fcrnandes et al. 2007) with predications from 
population synthesis models. 

For 21 galaxies in our sample, their stellar masses have 
been derived in ZGG11 by using optical spectra and stel- 
lar population synthesis code, STARLIGHT 5 (Cid Fer- 
nandes et al. 2005, 2007). During the fitting process, 
the Bruzual & Chariot (2003) models were adopted, and 
the spectra of the simple stellar population library were 
computed with the Salpcter (1955) initial mass function 
(IMF), Padova 1994 models, and the STELIB library (Le 
Borgne et al. 2003), by limiting the metallicities to 0.4Z Q 
for this metal poor sample. Stellar masses obtained with 
this method are represented by 2 (using integrated spec- 
tra; Moustakas & Kcnnicutt 2006) or 3 (using SDSS fiber 
spectra, and the fraction of the total galaxy luminosity 
in the SDSS i band inside the fiber had been used to 
correct for the aperture effect) in column 10 of Table 1. 
Please refer to ZGG11 for more details. 

For the rest galaxies, we use the M+/L ratio-color rela- 
tions given by Bell & de Jong (2001), which are derived 
by using the Bruzual & Chariot (2003) model with a 
Salpeter IMF, to convert luminosity to stellar mass. Due 

5 STARLIGHT & SEAGal: http://www.starlight.ufsc.br/. 



to the ongoing starburst activities in BCDs, we prefer to 
use the near infrared luminosity (Lk s ) and M+/Lk, to 
calculate the stellar mass. This is because that the K s - 
band luminosity is relatively less sensitive to recent star 
formation activity and dust extinction, and can well ap- 
proximate the stellar mass of a galaxy. The M+/LK a 
ratio-color relation, as given by Bell & de Jong (2001), 
is applied to correct for the dependence of M+/Lk s on 
star formation history. For our metal-poor BCDs, we 
adopt the result derived with metallicity Z = QAZq. 
In addition, 15 galaxies in our sample have both optical 
spectra and NIR photometric data. Therefore, we can 
compare the derived stellar masses. We find that the 
NIR-based stellar masses are systematically larger (0.11 
dex with a dispersion of 0.23 dcx) than spectrum-based 
ones. Therefore, we subtract this systematic difference 
from the NIR-based stellar masses. The results obtained 
with this method are also listed in Table 1, as represented 
by 1 in column 10. 

Eleven galaxies of our sample have neither optical spec- 
tra nor NIR photometry. For these galaxies we use R- 
band data and the M^/L^-color relation in Bell & de 
Jong (2001) to calculate their stellar masses. The i?-band 
data might be more affected by recent star formation and 
extinction than NIR data. To check to what extent, we 
compare the i?-band-based stellar masses with spectrum- 
based ones for 21 galaxies. We find that, in general, R- 
band-based stellar masses are 0.16 dex (with a dispersion 
of 0.27 dex) larger than spectrum-based ones. Hence, a 
correction of -0.16 dex is applied to the i?-band-bascd 
stellar masses, as represented by 4 in column 10 of Table 
1. 

As shown above, there exist some offsets between the 
stellar masses derived with different methods. The most 
possible reason causing this offset is the star formation 
histories (SFHs), as discussed in Gallazzi & Bell (2009). 
For galaxies with busty (or continuous) SFHs, such as 
BCDs, the one color-based mass-to-light ratio is gener- 
ally overestimated by 0.1-0.2 dex (Gallazzi & Bell 2009), 
which agrees well with our results. The dust might also 
play a role in this discrepancy since we do not correct 
for the internal extinction for the optical data. How- 
ever, dust extinction decreases the luminosity and red- 
dens the optical color. To first-order, these effects cancel 
out, and the estimated stellar mass of a galaxy is al- 
most unaffected even if a moderate internal extinction 
of E(B — V) = 0.2 mag is assumed (using the Cardelli 
et al. (1989) extinction law). The uncertainty in the de- 
rived stellar mass depends majorly on the uncertainties of 
models, galaxy evolution prescription, and SFHs, which 
are all of order of 0.1-0.2 dex or less (Bell & de Jong 
2001). Therefore, we adopt a statistical uncertainty of 
0.3 dex for our results. 

The stellar masses for the dl sample are derived with a 
two-component ( "young" plus "old" ) method (Lee et al. 
2003). To check whether this method is consistent with 
ours, we calculated the M+/Lb using the (B — V) color 
and M^/Ls-color relation given by Bell & de Jong (2001) 
for the dl sample, and then compared them with those 
derived with the two-component method. We found that 
the one-color based M*/Lb is systematically larger by 
0.14 dex, with a scatter of 0.11 dex, than that derived 
with the two-component method. This is generally in ac- 
cordance with the result described above. Therefore, we 
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conclude that the methods used to determine the stellar 
masses for the dl and BCD sample are consistent. 

2.3. Gas Mass and Fraction 

The total mass of atomic hydrogen, Mhi, in solar 
masses, can be derived from the Hi 21 cm observations 
by using the following equation (Roberts 1975; Roberts 
& Haynes 1994): 

il/ H i = (2.356 x lO 5 )^!^ 2 , (1) 

where F21 is the 21 cm flux integral in Jy km s _1 and D 
is the distance in Mpc. To account for helium and other 
metals, the total gas mass in solar masses is given by 

M gas = Mhi/X (2) 

where X is the fraction of the gas mass in the form of 
hydrogen, here adopted to be 0.733 (Lee et al. 2003) 
in consideration of the similar range of metallicity and 
the consistency between the dl and BCD samples. In 
Table 1 we include the logarithm of the total gas mass 
for our BCD sample, calculated from F21 (collected from 
literature; see references in Table 1) and the distance 
modulus included in the second column. 

The gas fraction, which is a natural parameter for 
quantifying to what degree a galaxy has exhausted its 
available fuel supply, is defined as 

_ Mgas _ M gas 

Af to t M gas + M* 

where the total mass in baryons, Af tot , is taken to be the 
mass of gas and stars. The gas fraction is a distance- 
independent quantity, because the gas mass and stellar 
mass are both derived from electromagnetic fluxes. De- 
rived gas fractions for the sample of BCDs are summa- 
rized in Table 1. For the BCDs in our sample, about 75% 
have their /1 > 0.3, and the median value of fi is 0.57. 

3. RESULTS AND ANALYSIS 

3.1. Correlations Involving Stellar Mass 

3.1.1. Stellar Mass-Gas Mass Relation 

The total gas mass and stellar mass in our BCD sample 
both range from ~ 10 6 to ~ 10 10 M Q . Figure 2 shows 
the distribution of total gas mass as a function of stellar 
mass for E-type and I-type BCDs. The solid line is a 
reference line for the case when the two qualities are the 
same. We find that, almost all I-type BCDs lie above 
the 1 : 1 line, whereas more than half of the E-type 
BCDs just below this line, indicating that I-type BCDs 
tend to have larger gas masses than E-type BCDs at 
a given stellar mass. For I-type BCDs, except for the 
apparent merging system, Mrk 108, whose gas content 
is extraordinarily high, the gas mass well correlates with 
the stellar mass, namely galaxies with a larger stellar 
mass also contains higher gas masses. The E-type BCDs 
seem to have a similar M gas — M* relation when M* < 
10 9 M Q . While at the high-mass end (M* > 1O 9 M ), 
galaxies tend to have lower gas contents comparing to the 
lower mass counterparts. As a comparison, in Figure 2 
we also plot the data for field dls which are adopted from 
Lee el al. (2003). In this plot, field dls seem to distribute 
similarly to BCDs as a whole, but the M gas — M* relation 
of field dls has a smaller dispersion than that of BCDs. 
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Fig. 2. — Distribution of total gas mass for I- and E-type BCDs, 
and field dls, as a function of stellar mass. The data for dls arc 
adopted from Lee et al. (2003). Solid line indicates 1 : 1 ratio, 
and the cross in the bottom right gives the typical error bar of the 
BCD sample. 

Since we are interested in the potential for connec- 
tions/differences between different types of dwarfs, it is 
informative to consider the gas fraction of a galaxy. In 
Figure 3 we plot the number distribution of gas fraction, 
log /it, in E- and I-type BCDs, and field dls. As men- 
tioned above, I-type BCDs are generally more gas-rich 
than E-type BCDs. The separation between these types 
of BCDs is more striking in Figure 3. The log// for E- 
type BCDs shows a wide distribution and ranges from 
~ —1.6 to <~ 0.0, and has a median value of —0.37. How- 
ever, the distribution of log /1 for I-type BCDs is much 
narrower, and looks similar to that for field dls. I-type 
BCDs concentrate between -0.2 - 0.0 (14 out of 17), 
while relatively more field dls have — 0.6 ^ log// ^ —0.2. 
The median values of log /i are —0.12 and —0.19, for I- 
type BCDs and field dls, respectively. 

To quantify the significance of the difference or sim- 
ilarity of the gas content among these different types 
of galaxies, we perform a nonparametric statistical test, 
the Kolmogorov-Smirnov (K-S) test. The K-S test for 
the I-type sample against the E-type sample rejects the 
null hypothesis (that these two data sets can be derived 
from a same parent population) at a 99.95% significance 
level. On the other hand, the K-S tests give probabilities 
of 6.65%, 3.22% and 26.64%, for the dls vs I- and E- 
type BCDs, all BCDs as a whole, data sets, respectively. 
Therefore, we conclude that the gas fraction distribu- 
tions of I- and E-type BCDs are much different, whereas 
the gas fraction distributions of dls and BCDs seem to 
be indistinguishable. Nevertheless, we caution that such 
analysis is possibly limited by the small numbers of the 
samples. 

3.1.2. Stellar Mass-Gas Fraction Relation 

We show the gas fraction of different types of galaxies 
plotted against their stellar mass in Figure 4. At the low- 
mass end (M* < 6.3 x 10 8 M Q ), in general, the gas mass 
is roughly independent on the stellar mass (i.e. propor- 
tional to the baryonic mass). At higher masses, however, 
the relation exhibits a break and becomes considerably 
steeper since the gas has been more efficiently converted 
to stars in higher stellar mass galaxies. 
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Fig. 3. — Number distribution of gas fraction for different types 
of dwarf galaxies. I-type BCDs have a higher mean gas fraction 
than E-type ones and field dls. 

From Figure 4 we can see that most BCDs and field 
dls distributes similarly in the log /it-M* plot. However, 
at a given stellar mass, the gas fractions for E-type BCDs 
have a larger span and a lower mean value than those for 
Ftype BCDs and field dls, which confirms that E-type 
BCDs are generally more gas-poor. 

3.2. Correlations Involving Metallicities 

Star-forming dwarf galaxies are observed to be well 
mixed chemically throughout their high surface bright- 
ness components (Kobulnicky & Skillman 1996, 1997), 
having no significant gradients or inhomogeneitics in oxy- 
gen abundances measured from H n regions located at 
different galactocentric radii (e.g. Lee & Skillman 2004). 
In Figure 5 we show the metallicity distributions for E- 
and Ftype BCDs and dls. It seems that all of these dif- 
ferent types of dwarfs show a double-peaked distribution 
of metallicity. However, E-type BCDs have a larger me- 
dian oxygen abundance (~ 8.2) than I-type BCDs and 
dls (both ~ 8.0). In the following we explore the chemi- 
cal evolution for BCDs by investigating in detail the cor- 
relations between their oxygen abundance and gas mass, 
baryonic mass, and gas fraction. 

3.2.1. Metallicity versus Gas Mass 

Figure 6 shows the metallicity as a function of the total 
gas mass. Albeit with large scatters, a galaxy with a 
greater gas mass is generally more abundant in metals. 
If we look into different types of galaxies, the scatter will 
reduce much. A unweighted least-squares linear fit, using 
a geometrical mean functional relationship (Isobe et al. 
1990), to the two BCD subsamples, gives, 

12 + log(0/H) = (3.65 ± 0.84) 

12 + log(0/H) = (3.34 ± 1.29) 



(0.55 ± 1.42) logM ga 
(0.53 ±0.32) log M ga 



for E- and I-type BCDs respectively. The fitted trends 
are shown as (purple) solid and (blue) dashed lines in 
Figure 6, and the rms deviations of the data from the 
relationships are 0.35 and 0.34 dex in (O/H), for E- and 
I-type BCDs, respectively. While for the dl sample, we 
obtain the following (O/H) — M gas relation 

12 + log(0/H) = (4.53 ±0.36) + (0.40 ± 0.10) log M gas . 
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Fig. 4. — Distribution of gas fraction (log/i) as a function of 
stellar mass for BCDs and dls. The cross in the bottom right gives 
the typical error bar of the BCD sample. More massive systems 
have lower gas fractions. 
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Fig. 5. — Number distribution of metallicity for different types of 
dwarf galaxies. E-type BCDs have a higher mean metallicity than 
I-type ones and field dls. 



with a rms scatter of 0.26 dex in (O/H), and is shown 
with the (red) dot-dashed line in Figure 6. 

Although the slopes of the Af gas — Z relations for I- and 
E-type BCDs are almost the same, and there are larger 
scatters in these relations, it is clear that E-type BCDs 
have higher mean oxygen abundance than I-type ones at 
a given gas mass, suggesting that they are depleted in gas 
relative to I-type BCDs of comparable metallicity. Field 
dls have a shallower Mg as — Z relation than BCDs, but 
are mixed with BCDs in the plot. Combining this dif- 
ference with the fact that I-type BCDs are more gas-rich 
than E-type BCDs (see Section 3.2), it might indicate a 
possible difference in potential for the progenitors and/or 
evolutionary phases of E- and I-type BCDs. 

3.2.2. Metallicity versus Baryonic Mass 

Figure 7 plots the oxygen abundance as a function of 
the total baryonic mass. For E-type BCDs, unlike Figure 
6, this figure shows a much better correlation between 
metallicity and baryonic mass. A geometric linear fit 
to the E-type sample (36 members) gives the following 
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Fig. 6. — Correlation between oxygen abundance and gas mass. 
The bottom right cross shows the typical error bar of the BCD 
sample. Albeit there is a large scatter, a galaxy with a greater 
gas mass is generally more abundant in metals. The slopes of this 
relation for I- and E-type BCDs are similar, and arc steeper than 
that of field dls. However, in general, E-typc BCDs have a higher 
mean metallicity than I- type ones at a given gas mass. 
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Fig. 7. — Correlation between oxygen abundance and the total 
baryon mass (Af gas + M*), with the cross showing the typical error 
bar of the BCD sample. 

relation 

12 + log(0/H) = (3.59 ±0.61) + (0.53 ± 0.11) log M tot . 

with a rms scatter of 0.25 dex in (O/H). The fitted rela- 
tion is plotted with a solid line. 

We also do the same fits to the I-type BCD and dl 
samples, which give 

12 + log(0/H) = (3.20 ± 1.22) + (0.53 ± 0.18) log M tot . 
and 

12 + log(0/H) = (4.63 ± 0.33) + (0.38 ± 0.06) logM tot . 

with rms scatters of 0.31 dex and 0.22 dex in (O/H), for 
the I-type BCDs and dls, respectively. The dashed and 
dot-dashed lines superposed in Figure 7 are the respec- 
tive fitting results for I-type BCDs and dls. 

These results show that the baryonic mass-mctallicity 
relations for I- and E-type BCDs have similar slopes, 
which are a bit steeper than that of field dls. However, 
E-typc BCDs appear to lie at a higher oxygen abundance 



comparing with I-type ones at a given baryonic mass, 
which again seems to suggest that E- and I-type BCDs 
might be at different evolutionary phases and/or have 
different progenitors. 

In fact, BCDs also have steeper L-Z and mass- 
metallicity (M-Z) relations than dls (ZGG10; Zahid ct 
al. 2011). Kunth & Ostlin (2000) show that there exist 
offsets between the B band L-Z relations of BCDs, dls 
and dEs/dSphs, in the sense that at a given B luminosity 
it becomes relatively more metal-deficient following the 
dEs/dSphs, dls and BCDs sequence. This offset persists 
even only comparing the metallicities of the old popula- 
tions in dls and dEs/dSphs (Grebel ct al. 2003). Woo 
ct al. (2008) find the same phenomenon in the M-Z re- 
lations of dls and dEs/dSphs. In the current work, how- 
ever, we find that dls distribute in between the two types 
of BCDs in the baryonic-metallicity relation. Therefore, 
the evolutionary connections, if exist, between dEs, dls, 
and BCDs may be complicated. 

3.2.3. Metallicity versus Gas Fraction 

The relation between metallicity and gas fraction can 
be used to check different chemical evolution models. 
The closed-box model predicts a linear correlation be- 
tween the logarithm of the inverse gas fraction \j \x and 
the oxygen abundance by number, n(0)/n(H), namely 
(e.g. Lee et al. 2003) 

12 + log(0/H) = 11.3 + logyo + loglog(l/ M ) (4) 

where yo is the yield of oxygen and we have assumed that 
the fraction of gas in the form of hydrogen is X = 0.733 
(solar value). 

A comparison of the observed oxygen abundances to 
those predicted by closed-box models (thick solid line) is 
shown in Figure 8, which indicates that few of the galax- 
ies in these BCD and dl samples agree with a closed-box 
chemical evolution. The vast majority of galaxies ap- 
pear to have oxygen abundances less than that expected 
by closed-box models (with a Salpeter IMF), which im- 
ply varying IMF, and/or gas inflows/outflows in these 
galaxies. It is also interesting to find that I-type BCDs 
and dls occupy a similar region in this plot, whereas E- 
type BCDs extend to the most metal-rich and gas-poor 
region. This indicates that (1) these galaxies follow a 
similar chemical evolution path (e.g. experiencing simi- 
lar gas flows), but at different phase (e.g. E-type BCDs 
are more evolved than dls and I-type BCDs); or (2) E- 
type BCDs experience different gas flows from dls and 
I-type BCDs (see below). 

The effects of inflows and/or outflows have been de- 
scribed by many authors (e.g., Edmunds 1990; Dalcan- 
ton 2007; Erb 2008; Matteucci 2008; Spitoni et al. 2010). 
In particular, Dalcanton (2007) shows that inflows are 
very efficient in reducing the effective yield, y e ff, for 
gas-poor galaxies while outflows are more effective in 
changing y e ff in gas-rich galaxies. Erb (2008) and Mat- 
teucci (2008) introduce a model in which galaxies have 
inflows/outflows proportional to the star formation rate 
(SFR). This model adopts the instance recycling and 
mixing approximations, and the infalling gas is assumed 
to have no metal, while the enriched outflows is con- 
sidered to have the same metallicity as the interstellar 
medium (ISM). In this model, we can change two free 
parameters, the amount of outflow f a and inflow fa in 
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Fig. 8. — Observed oxygen abundance plotted as a function of in- 
verse gas fraction for E- and I-type BCDs and field dls. The typical 
error bar of the BCD sample is indicated by the bottom right cross. 
Superposed lines are expectations from different models of chemical 
evolution. The thick solid line shows the result from a closed-box 
model with an oxygen yield approximate for low-metallicity galax- 
ies with a Salpeter IMF (y = 0.0074; Meynet & Maeder 2002). 
The thin solid (dashed) lines are pure outflow (inflow) models, with 
an effective yield 1 (black) and i (blue) of the true yield. The rest 
(noncontinuous) lines show the results of applying the model in 
Frb (2008) with values of f , fi and yo in the labels. The labeled 
galaxy, Mrk 108, which is an apparent merging system, appears to 
be interpreted by none of these models. 

unit of the SFR of the galaxy, to reproduce the observed 
oxygen abundance. Two more parameters can be var- 
ied, if necessary, i.e., the true yield y (i.e. IMF) and the 
fraction a of mass still locked in stars. 

In Figure 8 we also plot the predicted metallicity-gas 
fraction relations from this model using different param- 
eters. We adopt yo as 1 and | of the true oxygen 
yield predicted by the closed-box model with a Salpeter 
IMF( yo ,truc = 0.0074; Meynet & Maeder 2002), and 
vary fo an d fi to reproduce the dependence of oxygen 
abundance on the inverse gas fraction. The parameter 
a is a slowly decreasing function of time, as more stars 
leave the main sequence, but can vary in quite a small 
range. It is expected to be ~ 0.76 at t = 10 9 yr and 
~ 0.69 at t = 10 10 yr (Bruzual & Chariot 2003) for sub- 
solar metallicity galaxies with a Salpeter IMF. Given the 
age of these galaxies (see, e.g. ZGG11), we have used 
a = 0.7, in despite of that the results are not very sensi- 
tive to the chosen value. 

For a Salpeter IMF, as shown in Figure 8, the pure in- 
flow model (black dotted lines) with a variable rate can 
not explain the observed metallicity-gas fraction relation, 
(especially for E-type BCDs), which is consistent with 
the finding in Spitoni et al. (2010). Small inflow rate pre- 
dicts a higher metallicity than the observed one, whereas 
high inflow rate can not reach the gas-poor region. Pure 
outflow model (thin solid lines) or the combination of 
outflow and inflow models (dashed and dot-dashed lines) 
can reproduce the observed distribution, but the pure 
outflow model may need very high flow rate. The best 
solution could be a combination of gas flows and a vari- 
able IMF. 

4. DISCUSSION 

In this section we will mainly discuss the possi- 
ble evolutionary connections between BCDs, dEs and 



dls. There have been several evolutionary scenarios, 
both from observational and theoretical aspects, between 
BCDs, dEs, and dls (see Papaderos et al. 1996 for a re- 
view). However, the question of whether there is an evo- 
lutionary connection between these three types of dwarf 
galaxies have not been achieved a solid answer yet. As 
mentioned in Section 1, some apparently contrary results 
are present in previous studies on the basis of observa- 
tional measurements. 

For the BCDs in our sample, metallicity correlates with 
gas mass and total baryonic mass, in the sense that more 
massive systems are more abundant in metals. For these 
two correlations, I-type BCDs show similar slopes to E- 
type ones. However, I-type BCDs are generally more 
metal-poor and gas-rich than E-type ones at a given 
gas/baryonic mass. Although field dls seem to have a 
shallower slope relative to BCDs, they are well mixed 
with BCDs in the distributions of M gas — Z and M to t ~ Z 
plots. Moreover, in the plot of metallicity-gas fraction re- 
lation, I-type BCDs generally occupy a region differing 
from E-type ones. 

Our results may shed some light on the evolutionary 
connection between the three types of dwarf galaxies. 
I- and E-type BCDs might be at different evolutionary 
phases and/or have different progenitors. For example, 
if some E-type BCDs are formed from dEs (e.g. Silk 
et al. 1987), whereas I-type ones are formed from dls, 
then this part of E-type BCDs should be more evolved 
than I-type ones. This scenario is also supported by the 
fact that many E-type BCDs show structural properties 
compatible with dEs (Papaderos et al. 1996; Gil de Paz 
& Madorc 2005). In addition, field dls show chemical and 
gaseous properties in between I- and E-type BCDs, which 
suggests that they could be connected to both types of 
BCDs (e.g. Thuan 1985; Davies & Phillipps 1988). 

As discussed in Papaderos et al. (1996), the in- 
terchanges between dls, BCDs and dEs need gas out- 
flows/inflows to get involved, which is very consistent 
with our result that most BCDs should have experienced 
gas flows during their lives. Indeed, observations have 
provided (direct) evidence for gas flows in BCDs (e.g., 
Papaderos et al. 1994; Martin 1996, 1999; Martin et al. 
2001; Elson et al. 2009; van Eymeren et al. 2009; Lopez- 
Sachez et al. 2010). Using X-ray data, Papaderos et al. 
(1994) estimate a galactic wind flowing at a velocity of 
> 1000 km s" 1 for the BCD VII Zw 403, whereas the es- 
timated escape velocity for this galaxy is only 50 km s _1 . 
Martin et al. (2001) have shown that the galactic wind in 
NGC 1569 carries away almost all the oxygen produced in 
the current starburst. More recently, Elson et al. (2009) 
present their Hi observations for NGC 2915, and find 
a gas component clearly lagging in velocity relative to 
the main gas disk of the galaxy. This might be evidence 
for gas accretion from inter-galactic space onto the outer 
galaxy disk. Van Eymeren et al. (2009) report their de- 
tailed kinematic studies on the neutral (H i) and ionized 
(Ha) gas components in a nearby BCD, NGC 4861. The 
authors find that both in Hi and Ha the galaxy shows 
several outflows. Using multi-wavelength data, another 
work by this group (Lopez-Sachez et al. 2010) have stud- 
ied the environment of a sample of nearby BCDs, and 
shown that all BCDs present evident interacting features 
in their H I component despite the environment in which 
they reside. In particular, IC 4870, which has an ellipti- 
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cal outer component, is accreting H I clouds. 

5. SUMMARY 

We have presented a detailed studies on the gaseous 
and chemical properties of a sample of 53 BCDs, which 
is selected from G03/ZGG10 sample. According to their 
morphological information, this sample of BCDs is di- 
vided into two subsamples: E- and I-type BCD samples. 
An E-typc BCD contains an elliptical outer component 
while an I-type BCD has an irregular outer halo. 

For both types of BCDs, various correlations among 
metallicity, gaseous, stellar and baryonic mass, and gas 
fraction, are investigated. We find that in general, (1) 
more massive systems are more metal-rich, (2) a galaxy 
having a larger stellar mass also contains more gas, but a 
smaller gas fraction, (3) I- and E-type BCDs have similar 
slopes of the M gas — Z and M to t — Z relations, but I-type 
BCDs are generally more gas-rich and metal-poor. Based 
on the comparisons of these correlations, we suggest that 



I- and (at least a part of) E-typc BCDs might be likely 
at different evolutionary stages and/or having progeni- 
tors. By comparing the metallicity-gas fraction relation 
for our sample with those predicted by chemical evolu- 
tion models, we conclude that most BCDs should have 
experienced some gas inflows/outflows and/or mergers. 
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